Laser cooling all the way down to molecular condensate 
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Numerical simulations show that laser cooling of fermions 
on the repulsive side of the Feshbach resonance can sympa- 
thetically cool molecules well below their condensation tem- 
perature. 

PACS: 03.75.Ss, 39.25.+k, 03.75-b 

Fermi superfluids are in the center of interest of 
recent studies in physics of ultracold quantum gases. 
They provide perhaps the most promising model sys- 
tem to study the superfluidity, and the Bardeen-Cooper- 
Schrieffer (BCS) theory in the limit of strong interac- 
tions, when the size of Cooper pairs becomes comparable 
or even smaller than average interparticle distance, and 
the mean field description breaks down [1] . In cold atom 
systems the use of Feshbach resonances permits to tune 
the s-wave scattering length between the atoms with op- 
posite "spin" [2], assuring in this way a high flexibility 
of the systems. This technique allows in particular to 
drive the system through the BCS-BEC crossover when 
the scattering length diverges and changes sign from neg- 
ative (attractive) to positive (repulsive). During this pas- 
sage, the state of the system continuously evolves from 
the BCS superfluid of weakly correlated Cooper pairs, 
through the pseudogap regime, and then the unitarity 
limit, when the scattering length diverges, to the Bose- 
Einstein condensate (BEC) of diatomic molecules. By 
now several successful experiments with molecular con- 
densates were performed [3] which opened a way to first 
successful experiments with Fermi superfluids [4] . 

Earlier theoretical studies demonstrated that laser 
cooling allows to reach the degeneracy in Fermi systems 
[5], and to go all the way down to the Fermi superfluid 
[6] on the attractive side of the Feshbach resonance. The 
purpose of this Communication is to show that laser cool- 
ing offers also an alternative technique to reach molecu- 
lar condensate on the repulsive side of the resonance. 
Moreover, it allows to post-cool the condensate prepared 
via evaporative cooling down to even lower temperatures. 
We have developed here a fully self-consistent quantum 
kinetic theory of molecular condensate creation via laser 
cooling 1 . The results are very promising: despite the fact 
that condensate formation shifts the energies of quasi- 
particle excitations with respect to the equally spaced 
harmonic oscillator levels, and at the same time reduces 



1 similar to self-consistent theories of the BEC growth [7] 



fermionic population exposed to laser cooling, the cool- 
ing remains efficient well below the condensation temper- 
ature of 0.182V down to at least Q.OQTp. Optimization of 
the cooling protocol could probably significantly reduce 
this temperature. 

We consider system containing two species of fermions 
with "spin up and down" , molecules which are in their 
bound states, and fermions in an excited state. We work 
on the repulsive side of the Feshbach resonance, where 
interactions strengths between spin up/down fermions 
and molecules, as well as the binding energy of the 
molecules are simple universal functions of the effective 
s-wave scattering length a between spin up and down 
fermions [8]. Close to the resonance the binding en- 
ergy of molecules v = — 2-* can be made less than 
the condensation temperature ksTc for the molecules. 
As a result at all temperatures greater than Tc popu- 
lations of molecules and fermions are comparable. Be- 
low Tc, a fraction of molecules condenses leaving room 
in molecular phase space for more fermions to bind into 
molecules. Fermionic population is then shrinking, but 
it remains still comparable to the shrinking population 
of non-condensed molecules all the way down to the dis- 
sociation temperature, Tjj = f/feg- Below Try, fermionic 
population is suppressed, but at the same time almost all 
molecules are already condensed. Thus, close to the Fes- 
hbach resonance, when Tr> < Tc, there are favorable con- 
ditions for efficient sympathetic cooling of molecules by 
collisions with spin up/down fermions at all temperatures 
of interest. The fermions in turn are subject to laser cool- 
ing, and provide the heat reservoir, which is comparable 
to the cooled system (non-condensed molecules) . Popula- 
tions of the spin up and down fermions are the same. Mu- 
tual s-wave interactions between these two populations 
and the population of molecules lead to thermalization on 
a time scale, which is faster than the rate of laser cool- 
ing. At the same time three-body collisions keep the sys- 
tem close to chemical equilibrium between molecules and 
fermions. The system remains thus in a state of quasi- 
equilibrium with slowly decreasing temperature. This 
thermal state is described by the generalized BCS the- 
ory or boson-fermion model [9] which includes BEC of 
molecules described by the Gross-Pitaevskii equation. 
The model takes into account coherent tunneling be- 
tween molecules and pairs of spin up and down fermions. 
Thermal excitations of the molecular condensate are de- 
scribed by the bosonic Bogoliubov theory. Structure of 
the ground state of fermions and its fermionic excitations 
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follows from the set of Bogoliubov-de Gennes equations 
together with self-consistent definitions of the gap func- 
tion and the mean field potential. On the repulsive side 
of the Feshbach resonance the leading contribution to 
the gap function comes from the molecular BEC through 
the coherent coupling between molecules and pairs of 
fcrmions. In particular the gap function experienced by 
fcrmions is proportional to the condensate wave func- 
tion. Laser cooling of fermions was described in detail 
in Ref. [5]. Coherent laser excites atoms from, say, the 
spin down ground state to the excited state, and sponta- 
neous emission brings them back to the spin down ground 
state. Frequencies of the laser assure that the generalized 
Raman cooling takes place [10]. In this Communication 
we take into account interactions between fermions and 
describe the cooling process in terms of instantaneous 
fcrmionic and bosonic (molecular) Bogoliubov quasipar- 
ticles, whose eigenfunctions and eigenenergies are self- 
consistcntly updated during the evolution. We work in 
the "festina lente" limit to avoid reabsorbtion effects [11], 
and also employ spherical symmetry and ergodic approx- 
imations [6,7]. 

The boson-fermion model [9] is defined by the Hamil- 
tonian 



H= Jd 3 r 



a=+, 



tyLi>+i>- + h - c - + 

gi>\.ft_i>-i>+ + 0.3 g^\J)\J) m i} m + 

where the fields tp + and describe fermions with spin 
up and down, respectively, "tjj e corresponds to fermions 
in the excited state, and ip m is the bosonic molecular 
field. Here TL a = -^-V 2 + V a (r) - fi is a single particle 
Hamiltonian, m a is atomic mass m or molecular mass 
2m, V a (r) = ^m a Q 2 r 2 is the harmonic trap potential, \i 

is chemical potential, g — 47Th a is fermion-fermion in- 
teraction strength with a large positive s-wave scattering 
length a, A is a coupling between pairs of atoms and 

molecules, and v = r is the molecular binding cn- 

ergy. Close to the Feshbach resonance, when g is large, 
we can neglect interactions with excited atoms. 

In the mean field approximation, which closely follows 
the BCS theory, the quartic fermion-fermion interaction 
term in the Hamiltonian (1) is made quadratic by re- 
placing products of operators by their averages in all 
possible ways. The averages are the mean field poten- 
tial W{r) = 1.2 g Ul{r)^ m (r)\ +g (^(f)^ ± (r)\, and 



the anomalous pairing potential P(r) = \ij+{r)t}j-{r) 
which is mixing fermions with spin up and down. In 
the mean-field approximation the molecular field tp m is 



replaced by a c-number condensate amplitude <j> that ful- 
fills the stationary Gross-Pitaevskii equation 



(2 M - v)<f> = H m <t> + 0.6 <?|0| 2 c/> + XP . 



(2) 



Here (f> is normalized to the number of condensed 
molecules. The mean-field equations for fermions are 



lh dt 



H± ip± + W(r) i>± =F A(r) $ 



(3) 



with an effective gap function A(r) = gP{r) + \<j>(r) . 
These equations are mixing spin up and down compo- 
nents but they can be "diagonalized" by the Bogoliubov 
transformation t/j±(r) = £ m &m,± u m (f) T &L,=f "m( r 1. 
with fermionic quasiparticle annihilation operators 6 m ,±. 
The quasiparticle modes (u m ,v m ) fulfill the Bogoliubov- 
de Gennes equations with positive energies ui m , 



(dm U m 



H± U m + W Um - A V m , 

H± v m - W v m - A* u m . (4) 



In a thermal state with inverse temperature (3 the av- 
erage occupation numbers of quasiparticle states are 
N m = (bin ,±bm,±) = [exp(/3w TO ) + l]" 1 . Equations (4) 
are solved together with the self-consistency conditions 



A(f) = A0(f) + gJ2 [1 - 2iV m ] u m (r)v* m (r) 



W(f) = 1.2 g(Tpl(r)i, m (r)) + 



{1 - N m )\v m {r)\ 2 + N m \u m (r)\ 



(5) 



(6) 



by successive iterations with the chemical potential \i 
adjusted to keep the total number of atoms constant. 
The average molecular density {4> m (^)i ; m(r)} is a sum 
of the condensate density |</>(r)| 2 plus density of non- 
condensed molecules. After every iteration the Gross- 
Pitaevskii equation (2) is solved by relaxation. The to- 
tal number of atoms includes atoms bound into non- 
condensed molecules depleted from the condensate by 
thermal and quantum fluctuations. We estimate the 
number of non-condensed molecules using the approx- 
imate bosonic Bogoliubov modes in the Thomas-Fermi 
limit [12]. The ultraviolet divergence in Eq.(5) is regu- 
larized by the quickly convergent method of Ref. [13]. 

Excitation of atoms from the spin down state to the 
excited state is described by the Hamiltonian 



-I 



d 3 r 



i^e 4 ^> e t (r)V'_(r) + h.c. - <5^]4 



(7) 



driving coherent oscillations with Rabi frequency SI and 
laser detuning 6. The excitation is accompanied by the 
spontaneous emission described by a superoperator 
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£ P = 7 E U mk^\n-^k]p + V mk V[b mi+ e k ]p (8) 

with a spontaneous emission rate 7. Here, the Lind- 
blad superoperator is T>[A]p = ApA^ — \A^ A — 
\A^Ap, and the matrix elements are e.g. U mk — 
J dOfc W(f2fe)|u m fe(fc)| 2 with the spontaneous emission 
pattern W(fifc) and the generalized Frank-Condon fac- 
tors u mk (k) = J d 3 r e lkr wl(r)u m (r) . Here, Wfe(r) is the 
fc-th eigenstate of the harmonic oscillator and e k is an 
annihilation operator of an excited atom in this state. 

Adiabatic elimination of the excited state (cf. [14]) 
leads to the kinetic equations for the occupation num- 
bers taking into account the effects of Fermi statistics 

^ = E (1 - N ^-) N n,~ " (» ~ »») + 

n 

C mn (1 - JV m ,_) (1 - JV n ,+) - A nTO 7V mi _AT„ ;+ , 

n 

C„ m (1 - A m ,+) (1 - N ni -) - A mn N m .+N n ^ . 

They describe relaxation of ± quasiparticles (r^i_ n ), and 
creation/annihilation of pairs of + and — quasiparticles 
(C nm and A nm ). The transition rates are e.g. 



o 2 

r(-) = — V — 



l 2 Urnk\Unk(k L )\ 2 



(9) 



Here, w£ is the energy of the fc-th harmonic oscillator 
state and 7& is approximate spontaneous decay rate of an 
excited atom in this state, j k = 7 J2 m U m k(l - N m _) + 
V mk N mt+ , see Appendix A. 

Laser cooling drives average occupation numbers N m ,± 
of fermionic quasiparticles out of thermal equilibrium. 
At the same time interactions between fermionic and 
bosonic quasiparticles drive the system towards ther- 
mal equilibrium. Numerical simulations in Ref. [5] show 
that thermal relaxation remains efficient all the way 
down to the temperature of 0.03 Tp. This justifies our 
assumption of fast equilibration to a quasi-equilibrium 
state. After a short period of laser cooling dt the oc- 
cupation numbers N mj ± go out of equilibrium where 
N m ,± = N m = \exp((3u) m ) + The initial total en- 

ergy E(T) = £ m ,± co m N m . ± + Z m "X, including 
bosonic Bogoliubov quasiparticles with frequencies iv^, 
and average occupation numbers — [exp(/3w^) — 
changes by dE = E m w m (iV m , + +JV m ,_-2iV m ). 
The relaxation brings the occupation numbers N m ^± to 
a new state of equilibrium at a temperature T + dT, 
but docs not change the total energy. The energy of 
the system at the new temperature E(T + dT) differs 
from the initial E(T) by dE T w 2^^ m w m^m(l - 
N m ) + ^ E m «) 2 N£(l + N%). Conservation of en- 
ergy in thermal relaxation means that dE — dEx- In our 



simulations we use this equality to find the new lower 
temperature T + dT, and then solve the Bogoliubov-de 
Gennes equations and the Gross-Pitaevskii equation self- 
consistently to adjust the Bogoliubov modes, the conden- 
sate wave function and the chemical potential to the new 
lower temperature. With the new Bogoliubov modes we 
calculate new transition rates, etc. 

Calculation of the transition rates is the most time 
consuming part of the numerical simulation. This is 
why we were forced to assume spherical symmetry. For 
spherically symmetric A(r), W(r) and <p(r) the Bogoli- 
ubov modes are ui n {r) Yi m (6,<p) and their energies uji n 
do not depend on m. The bosonic quasiparticle en- 
ergies in the Thomas- Fermi limit [12] become iof n = 
fl m (2n 2 + 2nl + 3n + I) 1 / 2 . To be consistent with the 
spherical symmetry we also made ergodic approximation 
that quasiparticle occupation numbers Ni n> ± do not de- 
pend on m. In other words, we assume fast equilibration 
within each quasiparticle energy shell which keeps the 
system in a spherically symmetric state. More on spher- 
ical symmetry in Appendix B. 

In our simulations we use 81 harmonic oscillator lev- 
els, and assume that there are N = 10660 atoms with 
spin up and down which in the non-interacting case at 
T = is equivalent to 39 filled energy levels. The fre- 
quency of the isotropic trap is Q a — 2ir 2400Hz, the same 
for all fermions and molecules. The scattering length for 
the interaction between the two species is a = 1200A, a 
realistic value for interactions between two spin states 
\F = 9/2,m F = 9/2) and \F = 9/2,m F = 7/2) of 
40 K near the Feshbach resonance. In the natural har- 
monic oscillator units this scattering length gives the in- 
teraction strength g — 4.7, and the molecular binding 
energy v = —7.3. The wavelength of the cooling laser 
is A = 720nm, and the laser detuning is 6 = —\2Ml a . 
In all simulations the frequency of Rabi oscillations is 
much less than the spontaneous emission rate, Q = O.I7, 
so that the average occupation of excited state remains 
small. In our simulations we neglect W(r). The main 
effect of the mean field potential is an overall shift of en- 
ergy levels which can be compensated by laser detuning 
5. The coherent tunneling rate in trap units is A = 1. 

In Fig. 1 we show temperature as a function of time 
for laser cooling with two values of 7. Similarly as on 
the BCS side of the Feshbach resonance the nonzero gap 
function A(r) induced by the molecular BEC does not 
affect much the efficiency of laser cooling. Moreover, the 
cooling remains efficient even below the dissociation tem- 
perature of 0.18Tp. This is not quite surprising because 
even at such low temperature there is a non-zero frac- 
tion of fermions induced by coherent tunneling from the 
molecular BEC, which can be cooled by the laser. The 
initial cooling is faster for 7 = 10 but the final tempera- 
ture of 0.06Tp is lower for 7 = 1.25. The best strategy 
to reach low temperatures in reasonable time is to do the 
cooling in a few stages starting with large 7 and finishing 
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with small 7, as shown in the inset in Fig.l. In Fig. 2 
we show the growth of the molecular condensate in the 
same simulations as in Fig.l. 

To summarize, our numerical simulations show that for 
reasonable experimental parameters consistent with the 
requirements of the Festina lente regime, it is possible 
to laser cool fermions on the positive side of the Fesh- 
bach resonance all the way down to the molecular BEC. 
Despite the pairing effects, which produce a gap in the 
quasiparticle energy spectrum and the heat capacity of 
non-condensed molecules, it is possible to reach at least 
0.06 Tp in a time of a few seconds. 
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FIG. 1. Temperature T/Tf as a function of time for 7 = 10 
(solid line) and 7 = 1.25 (dashed line). Here Tp = 40.5hQ a , 
and the observed Tc ~ 0.6Tf. The cooling is faster for 
7 = 10, but the lower final temperature of 0.06Tf is achieved 
for 7 = 1.25. The inset shows T/Tf in a two stage cooling 
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FIG. 2. The number of condensed molecules in the same 
simulations as in Fig. 1. 
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APPENDIX A: TRANSITION RATES 



Laser cooling is realized by the Hamiltonian 

H las = Ifi J d 3 r e in ^i>t(r)i>-(r) + h.c. 



(10) 



driving coherent oscillations with Rabi frequency ft between the excited state e and the state — , together with the 
spontaneous emission superoperator 

C = 7 J dip dcos9 W(9,if) J d 3 n d 3 r 2 e^ 1 "^ 

x 2^_(f 2 )tp e {r2)p 'ipt(f 1 )ip-(f 1 ) - ^i{f 1 )^-(f 1 )^_(f 2 )i> e (r2) p - P ^(ri)^_(fi)^i(r2)^e(r2) , (11) 

with an effective spontaneous emission rate of 2j and the fluorescence dipole pattern W(0, ip). Expansion in the 
eigenmodes wi (r) of the harmonic trap 



and the Bogoliubov transformation (see main text) give 



las 



ml 



u m i(k L ) 



h.c. 



(12) 



(13) 



Here the (generalized) Frank-Condon factors are 

u m l(k) = J d 3 r e lkp wUr)u m (r) , 
v* ml (k)= [d 3 re^ W ;(r)v; n (r). 



(14) 

In a similar way, and after rotating wave approximation for e-atoms and the Bogoliubov quasiparticles, the spontaneous 
emission term becomes 

£ P = "f^2 U ml (l.bln^il p e]b m - ~ ejej&m,-&m,- P ~ P ^l^m-Hn,-^ + 
ml 

7 E Vml ( 2 ^+ei p 4 b ln,+ - eleiKn,+t>m,+ P~P e\hV m ^b m , + ^j ■ (15) 

ml 

Here the spontaneous emission rates are 

Umi = J dipdcosO W{6,Lp) \u ml {k)\ 2 , 

Vmi = J dtpdcosO W{6,ip) \v ml (k)\ 2 . (16) 

Adiabatic elimination of the excited state e, similar as in Ref. [14], results in the following evolution equations for the 
occupation numbers 



dN, 



E r ™-« ( X - N ™.-) N ">~ - T n~2m (1 - ^m,- + C roB (1 - JV m ,_) (1 - iV„,+) - A nm N m ,_N n , + , (17) 

n 



dN m ,+ 



There are contributions from four different processes: 
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• relaxation of — quasiparticles 



r (-) =—Y^ l 2 Umi\u n i(k L )\ 2 

" 2 7 ^ \ S - U f + uj n + i 7 £ s U sl (l - iV s ,_ + 6 a , n ) + V sl N Si+ \ 2 ' 



relaxation of + quasiparticles 



r£L = %-Y. i 2 v nl \ Vml (k L W (20) 



'rn.s ) 



• creation of pairs of + and — quasiparticles 



mn 2 7 ^- \6-uf-u n + i^„ U sl (l-N s ^) + V sl (N s , + +S n , s )\ 2 ' 



annihilation of pairs of + and — quasiparticles 



A = ^!y- i 2 v m i\u nl (k L )\ 2 

mU *r 2 j'\5-u>f+u> n + iy'E a U B t(l-N Bt -+6 n , a ) + V Bl N Bt+ \ 2 ' 



APPENDIX B: SPHERICAL SYMMETRY 



Simulation of laser cooling with a non- vanishing pairing function P(f), and a Hartree potential gop(r) is a numer- 
ically hard problem. The main difficulty is that after every short period of laser cooling it is necessary to reiterate 
Bogoliubov-de Gennes equations in order to obtain a new self-consistent pairing function, Hartree potential, and the 
Bogoliubov modes (oj m , u m , v m ). With the new Bogoliubov modes, the transition rates A mn , C mn , T^ lt _ n have to be 
calculated again, and this is a very time-consuming operation. This numerical effort is reduced a lot with assumption 
of spherical symmetry. With spherically symmetric A(r) and p{r) the Bogoliubov modes can be decomposed into 
spherical and radial parts 

ui n (r) Y lm (e,cp) , (23) 
vi n (r)Y lm (6,<p) . (24) 

In a similar way harmonic oscillator modes become Wi n (r)Yi m (9, ip). Quasiparticle energies LUi n do not depend on m. 
We assume fast equilibration within each degenerate energy shell nl so that all occupation numbers in a given shell 
are the same iV|„,±. This assumption greatly simplifies the rate equations 

Chn 1 l 2 n 2 (1 _ ^lim ,-) (1 _ ^hn 2 ,+ ) ~ a l 2 n 2 hm ^Imi ,-^l 2 n 2 ,+ ! (25) 
^N knit+ = J2 hn2 lltl^l2n 2 (! - ^1-1,+) N hn 2 ,+ ~ T^l^ini i 1 ~ N hn 2 , + ) ^lim.t + 

ci 2 n 2 i ini (1 - ^iini.t) (1 _ Ni 2n2 _) - a hni i 2n2 Ni inii+ Ni 2n2 _ . (26) 
Here 7^, a, c are averaged transition rates, for example 

-y ( ~ } = 1 V r (_) (27^1 

il 1 n 1 ^l 2 n 2 _|_ \ / j lmim\<—l 2 n 2 m 2 " V ■ / 

ra\m 2 

This equation combined together with the definition of r( _ ) gives 

(-) _ 7^ \ - \ - Ui ( k L)\ 2 

ini ^ 2 " 2 2 ^ + !) J^ 2 J^n. V + 2n e - M) + + l 2 Rt nc ■ 



G 



Here the frequency of the intermediate excited state u)(l e + 2n e — fi) is measured with respect to the chemical potential 
H because the quasiparticle energies o>; 2 „ 2 are also defined with respect to /i. For a laser beam with a along the 
2-axis the indices m e and m 2 must be the same. Furthermore, it follows from the properties of Clebsch-Gordan 
coefficients that the sum J2 mi Ui in imil e n e m e does not depend on m e so that we can set e.g. m e = and 

-f£l 2 „ (Eroi Ul inimi l ene o) \ J2m 2 \ u l 2 n 2 m2l c n c m 2 

7 ' ini ^ 2 " 2 " 2(2/! + 1) fa {8-u;(l e + 2n e -n) + Lu n } 2 + 1 2Rl nF 



7 ft 



2 ( 2 *i + l) h l S - "(h + 2n e - M ) + c„] 2 + 7 2i?f 



inil e n e SUl 2 n 2 l eUe 



(29) 



The last form shows that the transition rate 7 ( ' can be constructed out of matrices SU, su and a vector i?. These 
elements can be expressed through even more elementary building blocks 



SU hniUne = ^2(jWu)l enJini S Ueh 



su l 2 n 2 l c n e — ^ ' 
m 

The more elementary building blocks are 



^2su e i 2m (jWu) Hene 



Sll.h = (21 + l)^tiy]((,(„ii>,0|i 1 ,m) 2 (i,i„0,0|i 1 ,0) i 

+ 1 



2L + 1 



sii e i 2 m = i\u + 1) y ^n^' /e ' °- m| ' 2 ' m) <*' °' 0|/2 ' 0) 

The line- width of an excited state l e , n e is approximately given by 

Rl e n e = X/^ 1 ~~ Ni n )SUl n l ene + Ni n SVi n i crli . 



(30) 



(31) 



The matrix SV is obtained from the matrix SU by a substitution u — > u. As the eigenfunction M; n (r) evolves in the 
process of laser cooling it is more efficient to express the matrix (jWu) through a static matrix (jWW) build out of 
harmonic oscillator modes 



(jWW) 



f . 

U e n e lmi I r 



dr ji{kr)Wi e „ e (r)Wi ini (r) 



(32) 



The static matrices s,S,jWW were prepared once, and stored then on a disk. They facilitate calculation of the 
temperature dependent matrices SU, su, R every time the wave functions u, v are updated. 
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